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A. PENGANTAR 

Transformasi Householder sering juga dikenal dengan istilah refleksi householder 
merupakan salah satu metode yang dapat digunakan dalam dekomposisi QR. Penamaan 
Householder sendiri berasal dari penemunya, yakni Alston Scott Householder. Dia seorang 
matematikawan Amerika yang ahli di bidang Matematika Biologi dan Analisis Numerik. 

Dalam tulisan ini saya akan memperlihatkan cara menghitung matriks Householder. 
Saya akan mulai dengan konsep tentang vektor yang perlu diingat kembali. 

Definisi 1 

Diberikan v, w e W , maka 

a. Panjang vektor v (disebut juga norma v) didefinisikan dengan persamaan 

IHI = V v i 2+v 2 2+ "" +v » 

b. Hasil kali titik (dot product) antara vektor v dan w didefinisikan dengan 

v .w = v ] w 1 +v 2 w 2 + --- + v„w„ 



Berdasarkan definisi di atas, maka hubungan dari norma vektor dengan hasil kali titik 
diperlihatkan oleh persamaan berikut 

v»v = ||v|| 2 atau ||v|| = yjv*v 
Umumnya, jika vef" maka v = (v 1 , v 2 ,- •, v„) . Namun dalam tulisan ini, saya akan 



memandang vektor v sebagai matriks kolom, yakni v 



mudah dapat dibuktikan bahwa v T v = llvll 



Akibatnya, dengan sangat 



B. PEMBAHASAN 

Definisi berikut ini akan memperlihatkan bentuk dari matriks Householder. Bentuk 
ini sebenarnya muncul dari konsep proyeksi orthogonal (tidak diberikan dalam tulisan ini). 

Definisi 2 

Transformasi Householder didefinisikan dengan persamaan 



v v 

Matriks H pada Persamaan (1) sering disebut matriks Householder. Kegunaan dari 
transformasi Householder adalah untuk merubah sebarang vektor tak-nol menjadi vektor 
yang semua elemennya bernilai nol kecuali untuk elemen pertama dari vektor tersebut. 



Dengan kata lain, untuk sebarang jeM" yang tak-nol maka 

"l' 





Hx = k 







= ke l , dimana k e R dan e x e W . 



(2) 



Dalam penulisan ini, saya akan mencoba menurunkan rumus untuk menghitung 
vektor v yang merupakan elemen penting dalam transformasi Householder. 
Langkah pertama, saya akan coba mengalikan matriks Householder dengan sebarang vektor x 
untuk melihat bentuk dari vektor v. Diperoleh 



f 



T \ 



Hx = 



1-2 



vv 

T 

v v 



X 



= x — 



= x — 



J 

f T \ 

2^ 

T 

V V V J 

f T \ 
T 

V v v J 



x 



= x-6v 

Gunakan hasil terakhir di atas dan Persamaan (2), sehingga diperoleh 
Hx = Hx => x-Ov-ke x 

=> 0v-x-ke, 



Selanjutnya, saya akan misalkan v = x- ae l . 
Hx = x- 
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^x-ae^ T x 
(x T -aefj. 



(x T -ae\^{x-ae^) 



(x-ae x ) 



[x-ae^ 



(x T x — ax 1 } 
xx — 2ax, + a 



2^x T x-ax l } 
x J x-2ax i + a 2 



Agar Persamaan (2) terpenuhi maka bagian 
Dengan kata lain, 
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2(x T x-ax^ 



x T x-2ax x +a 



2^x T x-ax^ 

x J x - 2ax, + a 2 
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harus sama dengan nol. 



2^x T x — ax 1 } 
x T x-2ax, +a 2 



= <=> 
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x T x-2ax, +a 
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Dari ruas kanan persamaan di atas, maka saya bisa menurunkan nilai dari a , sebagai berikut 
2^x T x — ax 1 ) 



= 1 => 2[x T x — ax 1 ) = x T x — 2ax x + 
xx- 2ax l +a v ' 

=> 2x T x-2ax l = x T x-2ax l +a 2 

T 2 

=^> x x = a 



a 2 



.2 



,2 



=^> a = ||x|| 
=> a = ±|x|| 

Jadi diperoleh v = x± \\x\\ e x . 

Masalah yang muncul adalah tanda mana yang harus dipilih? Apakah "+" atau "-" ? 
Untuk mengatasi masalah ini, saya akan andaikan bahwa x - x l e l , dimana x\ adalah elemen 

pertama dari vektor x. Jika x l > maka v = x — \x\e l akan bernilai nol, sebaliknya jika x x < 

maka v = x + ||x||e 1 yang akan bernilai nol. Oleh karena itu, persamaan v = x + ||x|| e x harus 
dirubah menjadi 

v = x + sign(x 1 )||x||e 1 (5) 
Perlu diketahui bahwa fungsi sign(xi) akan menghasilkan tanda("+" atau "-") dari x\. 
Namun jika x 1 = maka pilihlah tanda "+". 

Agar lebih jelas, berikut ini saya berikan contoh disertai hasil perhitungannya. 
Contoh 1 

^2 ' 

Tentukan matriks Householder untuk vektor x = 3 

-1 

Langkah pertama yang harus dilakukan adalah mencari vektor v dengan menggunakan 
Persamaa () di atas. Saya menggunakan aplikasi Matlab untuk menghitung dan diperoleh 
hasil sebagai berikut 
»x=[2;3;-l] 

x= » v = x + sign(x(l))*norm(x)*el 

2 

3 v = 

-1 5.7417 

3.0000 

»el=[l;0;0] -1.0000 

Dengan menggunakan vektor v ini maka bisa dihitung matriks Householder dan diperoleh 
hasil sebagai berikut 
» I=eye(3) 
1 = 

1 
1 o 
1 



» H=I-(2*v*v')/(v'*v) 
H = 

-0.5345 -0.8018 0.2673 



-0.8018 0.5811 0.1396 

0.2673 0.1396 0.9535 
Selanjutnya, saya akan coba mengalikan matriks H dengan vektor x untuk melihat hasilnya. 
» H*x 

ans = 
-3.7417 
-0.0000 
0.0000 

Dari hasil di atas terlihat bahwa semua elemen bernilai nol kecuali untuk elemen pertama. 

Contoh berikut memperlihatkan kasus x l = . 
Contoh 2 

Tentukan matriks Householder untuk vektor x 



Hitung matriks v dan diperoleh 
» x=[0;-3;4] 
x = 



-3 

4 

» v=x+norm(x)*el 
v = 

5 
-3 
4 

Selanjutnya hitung matriks H dan diperoleh 
» H=I-(2*v*v')/(v'*v) 
H = 

0.6000 -0.8000 

0. 6000 0.6400 0.4800 
-0.8000 0.4800 0.3600 

Jika saya kalikan matriks H dengan vektor x akan diperoleh 
» H*x 
ans = 

-5 





Catatan : 

1. Perintah norm(x) pada Matlab akan menghitung norma (panjang) dari vektor x. 

ii. Tulisan v kemudian diberi tanda kutip satu : v' memiliki arti v (vektor v ditransposkan) 

iii. Perhitungan standar pada Matlab menggunakan 4 digit di belakang tanda koma. 




C. PENUTUP 

Mohon maaf jika terdapat kekurangan ataupun kesalahan. Saran dan kritik dapat 
dikirim ke email saya : bernypebo@yahoo.co.id 



D. LAMPIRAN 



Dalam pembuktian di atas terdapat persamaan (vv r )x = (v r xjv . Saya akan buktikan untuk 
kasus v, x e E 3 , sebagai berikut : 



(vv T )x = 



3. 
V.V., 



[ V l V 2 V i] 
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"l "3 
V 2 V 3 



VjV 3 v 2 v 3 



X, 



v 2 l x l +{v l v 2 )x 2 +(v x v 3 )x 3 
(v l v 2 )x l +v 2 2 x 2 +(v 2 v 3 )x 3 
(v l v 3 )x l +(v 2 v 3 )x 2 +v 3 2 x 3 



(L. I) 



Di sisi lain, 
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v T x)v = 
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_ V 3 



= (v l x l + v 2 x 2 + v 3 x 3 ) 



v l 2 x l + (v l v 2 )x 2 + (v l v 3 )x 3 
(v l v 2 )jc l +v 2 2 x 2 +(v 2 v 3 )^ 

_(v l V3)x l +(v 2 V 3 )x 2 +V3 2 X3_ 

Berdasarkan (L.l) dan (L.2) maka terbukti (vv T )x = (v T xjv . 



(L.2) 



